Phase Separation under Ultra-Slow Cooling: Onset of Nucleation 



00 

o 
o 

in 



J. VollmeiB 

Fachbereich Physik, Philipps Universitat, 35032 Marburg, Germany 
Max Planck Institute for Dynamics and Self-Organization, Bunsenstr. 10, D-31013 Gottingen, Germany 

(Dated: March 15, 2008— 12:46pm) 

We discuss the interplay between a slow continuous drift of temperature, which induces continuous 
phase separation, and the non-linear diffusion term in the 4 -model for phase separation of a binary 
mixture. This leads to a bound for the stability of diffusive demixing. It is argued that these 
findings are not specific to the 4> 4 model, but that they always apply up to slight modifications of 
the bound. In practice stable diffusive demixing can only be achieved when special precautions are 
taken in experiments on real mixtures. Therefore, the recent observations on complex dynamical 
behavior in such systems should be considered as a new challenge for understanding generic features 
of phase-separating systems. 
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I. INTRODUCTION 

The kinetics of phase separation 
topic of continuous experimental @, S , 
theoretical interest fl3l 0, HH ■ Many of its characteris- 
tic features do not depend on specific materials, but are 
universal in the sense that they can be understood based 
on relatively simple model systems. For binary mixtures 
an appropriate model is the 4> 4 model 0, H [l6| . In the 
present paper we augment this model by terms taking 
into account a slow continuous drift of the temperature, 
in order to explore how the composition of coexisting 
phases keeps up with the driving due to the temperature 
evolution. 

Typically the demixing dynamics of binary fluids is 
idealized as a three step process: 

A. A rapid change of the temperature (pressure or other 
experimental control parameter) transfers the system 
from a single phase equilibrium into the region, where 
phase separation is expected [cooling in Fig. [Q. This 
is supposed to be sufficiently fast to neglect noticeable 
changes in the composition. 

B. The mixture decomposes into two coexisting macro- 
scopic phases [decomposition in Fig. [Q. In applica- 
tions this is supposed to be fast on the time scales of 
chan ging the temperature; in high precision experiments 
0) S Il7l [H, [lj| S3] a temperature jump is applied and 
the temperature is subsequentially fixed to carefully ob- 
serve the dynamics in this regime. 

C. Upon further slow cooling there is an exchange of 
material such that the composition of the coexisting 
phases follows the equilibrium ones up to a small time 
lag [smooth green line denoted segregation in Fig. Q]. 

A dimensional analysis of the diffusive relaxation time 
[2lL I22I S3 shows, however, that for experimentally 
relevant cooling rates the composition can not remain 
close to the equilibrium values unless temperature is 
ramped exceedingly slowly. Indeed, recent experimen- 
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FIG. 1: (Color online) Schematic phase diagram for the 
demixing of a binary mixture. The binodal is given in blue 
and the spinodal line in dark blue. A typical evolution of 
the composition (j> of the mixture upon constant cooling is 
indicated by the green line. The red line shows an oscillating 
demixing, as it may be observed when the diffusive decompo- 
sition becomes unstable. 
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ical studies [21|, |29( indicate that rather than slowly 
following the drift of the equilibrium composition, of- 
ten there is very complex dynamics observed in regime 
C. Among others secondary nucleation in large domains 
(double phase separation) [8|, oscillating bursts of mas- 
sive phase separation alternating with quiescent periods 
d dl [mjll EE ES E3, El and stationary convection 
patterns [2l[ have been reported. Except for the sta- 
tionary convection, all these scenarios give rise to an ex- 
plicitly time dependence of the composition like the one 
indicated by the red line in Fig. [TJ 

Despite the impressive experimental findings, the the- 
oretical understanding of such complex dynamical re- 
sponse to slow cooling is still at its infancy. To make 
a point, however, on the prevalence of the various com- 
plex patterns of phase separation the present paper ex- 
plores the critical cooling rate beyond which stable dif- 
fusive demixing can no longer be expected. 

The article is organized as follows. Section Sec. HTlre- 
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visits the relevant nonlinear diffusion equation describing 
the evolution of concentration in a slowly cooled binary 
mixture. In Sec. lIIII we consider the limit of a very narrow 
interface at the meniscus in order to analytically derive a 
phase portrait which gives us access to discuss stable and 
unstable solutions of the diffusion equation, and their bi- 
furcations. A generalization of this discussion to the case 
where the phases are not symmetric is given in App. [Al 
The subsequent Sec. |IV] introduces a mode expansion in 
order to gain insight in the time evolution of the profiles 
after bifurcations, and the role of higher order derivatives 
needed to properly describe interfaces of finite thickness. 
(Technical issues of its derivation have been delegated to 
App. [BP The main results are summarized in Sec. IVl 



II. NONLINEAR DIFFUSION EQUATION 

We describe the temporal evolution of the binary mix- 
ture by the free energy functional [l6[ 



F[<f>(x;t),t] = J dx 



F [4>(x;t)] + j (V0(x;i)) 2 



(1) 



where 0(x) characterizes the composition at position x, 
and the integration over x is performed over the volume 
of the system. The second term in square brackets repre- 
sents an energy penalty for steep changes of composition 
(macroscopically it gives rise to surface tension), and 



-Ml 



(2) 



is the free energy of the 4 -model, which describes the 
equilibrium phase behavior of the mixtures. In the 
present study we only deal with negative values of A, 
where there are two coexisting phases with composition 
±00 an d 0o = y—A/B. The temperature dependence 
of the equilibrium composition 0o(T) is due to the tem- 
perature dependence of the parameters A and B. 

The free energy ([2]) faithfully accounts for qualitative 
features of the demixing of fluid mixtures, even though 
more complicated expressions may arise for concrete sys- 
tems [4|, |l6j. The explicit time dependence in Eq. (TT]) 
accounts for the fact that we deal with a system sub- 
jected to a sustained change of temperature. Due to its 
temperature dependence the equilibrium composition 0o 
becomes explicitly time dependent in that case. 

In the absence of center of mass flow the interdiffusion 
of the components of the mixtures follows 

d t 0(x;i) = aV 2 /i[0(x;i)] 

where a is the thermodynamic coefficient for mass inter- 
diffusion, and the chemical potential /i[0] = <5(i*o[0] + 
K/2 [V0] 2 ) / 8<j) is the functional derivative of the free en- 
ergy density (TT]) with respect to 0. Consequently, the 
normalized composition ip = 0/0o evolves according to 



the nonlinear diffusion equation 
~D 



(3ip 2 - 1) V(p 



where 



'Jt(po 

00 



DfV>-^ (3) 



(4) 



accounts for the change of the equilibrium composition 
due to the changing temperature, S = (K/\A\y/ 2 is the 
width of a planar macroscopic interface between the co- 
existing phases fl(| . and D = —olA is the equilibrium 
diffusion coefficient {%. e., the one for a system with uni- 
form composition ip = ±1). Away from equilibrium the 
diffusion becomes nonlinear with a diffusion coefficient 
D (3ip 2 — l)/2 explicitly depending on composition. In 
the spinodal range — ip s = — 1/V3 < <p < l/v3 = tp s 
the diffusion coefficient even takes negative values, indi- 
cating the instability of the mixture against arbitrarily 
small fluctuations. 

In the following we focus on the simplest setting allow- 
ing us to discuss the breakdown of diffusion — general- 
izations are outlined in App. [A) 

(i) There are equal amounts of the components A and B 
of the mixture such that in equilibrium there is a macro- 
scopic A-rich phase coexisting with a B-rich one. 

(ii) The mixture is kept in a container of constant cross 
section and height 2A, and it is subjected to a gravi- 
tational field. As a consequence there is a macroscopic 
phase boundary at z = 0. The bottom and top walls of 
the container are at z = ±A, respectively. 

(iii) The flux through the bottom and top of the system 
vanishes. 

(iv) There is no flow in the horizontal directions such 
that d x ip = dytp = 0. 

(v) The temperature is changed in such a way that £ is 
constant (c/. [22|). 

(vi) The temperature and composition dependence of the 
equilibrium diffusion coefficient D may be neglected. 

Spacial scales are conveniently measured in units of A, 
and time in units of A 2 / D. The resulting dimensionless 
diffusion equation 



d t (p = d z 



3ip 2 - 1 



M 2 dfip-Nip (5) 



takes a universal form for all the mixtures, 
volves the dimensionless parameters 



N 
M 



£A 2 

S_ 

A 



It only in- 

(6a) 
(6b) 



The former characterized the strength of the driving by 
the ratio of the rate of generation £ of supersaturation 
over the rate of decay of supersaturation due to diffusion. 
The latter characterizes the importance of the boundary 
layer at the meniscus between the coexisting phases by 
the ratio of the interface width 6 over the system size A. 
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FIG. 2: (Color online) Stationary density profiles for M — and N = 0.1, 0.2, 0.3, 0.4, 0.50 (solid red lines from top to bottom), 
and the critical parameter N c = (v3— l) 2 — 0.5359 (dashed green curve): (a) stable solutions with <p'(z = 1) = 0; (b) unstable 
solutions, which all end at tp s — 3 _1//2 ~ 0.58. The dotted blue lines in part (a) show the approximation (|10p for N = 0.1 and 
N = 0.2, respectively. 



Equation §5§ ought to be solved subject to no-flux 
boundary conditions (iii) at the top and bottom of the 
system, viz. 







3^ 2 



1 



d z ip - M d 3 z ip 



(7) 



z=±l 



To ensure mass overall conservation we furthermore de- 
mand f{z) = —ip{z) [3l| . 

III. STATIONARY DIFFUSION PROFILES 

Except in the critical region (i e., the region very close 
to the maximum of the coexistence curve in Fig. [T]) the 
width 8 of the interface is always much smaller than the 
lateral extension A of a macroscopic system. Observ- 
ing that a temperature ramp rapidly moves the mixture 
out of this region, one may assume M = 6/ A -C 1 in 
systems of macroscopic dimension A. Consequently, the 
term with the fourth order derivative, which accounts for 
the effects of surface tension, is only relevant in situations 
where there are very rapid changes of composition. This 
only is the case close to z = at the interface between the 
coexisting phases, and for values of ip close to the thresh- 
old of spinodal decomposition (c/. below). Typically it 
has only negligible influence on the supersaturation pro- 
file ip(z;t). To gain insight in its dependence on N we 
therefore discuss the <p{z) profiles in the idealized case 
M — 0, where the diffusion equation 



d t ip = d z 



3<^ 2 



1 



d z tp 



Nip 



(8) 



will be subjected to the boundary conditions 

z = ± : p{z = 0) = ±1, (9a) 

Sip 2 -1 



z = ±l 



d 2 ip{z = ±l;t) = 0. (9b) 



The former requires that the the composition takes on 
the equilibrium value at the meniscus, and the latter that 



there is no flux through the top and bottom interfaces of 
the system. There are two possibilities to fulfill the latter 
boundary condition (Fig. [5]): 

(a) stable diffusion profiles are obtained for d z ip(z = 
±1) = 0; 

(b) when the composition at the top and/or bottom of 
the system exactly amounts to the spinodal composition 
tp(z = ±1) = ip s = 3 -1 / 2 the factor in front of the deriva- 
tive vanishes, and one obtains an unstable diffusion pro- 
file. 

These possibilities will be further explored in the remain- 
der of this section. The question in how far the find- 
ings based on this idealized setting still apply at a non- 
vanishing M will be addressed in Sec. IIVI 



A. Stable diffusion profiles 

For very weak driving, N 1, the composition ip 
never substantially differs from its equilibrium value, 
such that (3ip 2 — l)/2 ~ 1. In other words the diffu- 
sion coefficient hardly depends on ip, such that it may 
be approximated by its equilibrium value. In this case 
the boundary condition ([7]) at the top requires d z p = 0, 
and one easily checks that the stationary solution of the 
diffusion equation amounts to 



ip(z) = 



cosh 


Vn (z- 


-1)" 


cosh 


Vn 





(10) 



For N = 0.1 and N = 0.2 it is shown by a dotted blue 
line in Fig. Ela). Indeed the approximation appears to 
be very good as long as N < 0.1. 

For larger values of N the (^-profile deviates too 
strongly from unity to neglect the supersaturation de- 
pendence of the diffusion coefficient. In that case the 
solution of the diffusion equation (O can not be given 
in a closed form. Upon numerical integration one ob- 
serves, however, that close to z = it closer and closer 
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3- 




FIG. 3: (Color online) The values ip(z = 1/2) (thick green + 
connected by dahsed lines) and tp(z — 1) (thick red x con- 
nected by solid lines) of the stable (upper curve) and unsta- 
ble (lower curve) solution as a function of the dimensionless 
ramp rate N. The respective approximations by the two mode 
model Eq. (I18|l are indicated by corresponding small symbols. 

The corresponding small symbols show 



approaches the linear solution 



lfi(z) 



(11) 



Unless N = N c = (a/3 - l) 2 ^ 0.5359 this linear solution 
does not fulfill the no-flux boundary condition (|9b[) at 
the surface z — 0. Consequently, the solution eventually 
crosses over into an approximately parabolic profile with 
a minimum at z = 1 (solid red lines in Fig. 0Ja)). I n 
the limit N — > N c the stable solutions closer and closer 
approaches the linear solution 1 — ^/N c /3x, which for 
N = N c fulfills the boundary condition ([7]) because the 
ratio preceding d z ip vanishes. 



B. Unstable stationary profiles and a saddle-node 
bifurcation 



From a physical point of view the linear solution for 
N c has no flux across the top surface since the nonlinear 
diffusion coefficient (3ip 2 — l)D/2 vanishes at the spin- 
odal. There are several observation to be made about 
this solution: 

(i) Since the composition close to the top of the sys- 
tem approaches the value at the spinodal line, the linear 
profile is unstable against nucleation of domains of the 
minority phase close to top end. As a consequence the 
stable solutions of the diffusion equation approach an un- 
stable solution, when N — > N c . 

(ii) Besides being the limiting behavior of the stable so- 
lutions, the linear profile can also be viewed as a limit- 
ing case of a family of solutions that end on the spin- 
odal point at the top of the system (solid red lines in 
Fig. [Sfb)). All these solutions are unstable against nu- 
cleation of domains of the minority phase close to top 
end. 



FIG. 4: (Color online) Flow diagram of Eq. (S) with M = 
and N absorbed in the length scale (cf. text). The equilib- 
rium (ip = ±1) and spinodal (ip — ±l/v3) compositions are 
marked by vertical dashed green lines. Note that the lin- 
ear profiles (jlip run horizontally on lines (p'/y~N = ±l/v3, 
which are also indicated by dashed green lines. 



(iii) At the critical parameter value N c the family of sta- 
ble solutions of the diffusion equation collides with the 
unstable ones in a saddle node bifurcation. Subsequently, 
there is no stationary diffusive solution with only two 
domains in the vertical direction. In order to clearly 
make this point Fig. [3] shows the values <p(z — 1/2) and 
(p(z = 1) for the coexisting stable (upper curves) and 
unstable (lower curves) solutions. 

This result still holds for non- vanishing M, since dj<p 
vanishes for a linear profile. Moreover, it also does not 
rely on the specific dependence (3<p 2 — l)/2 of the (j) 4 
model. When the supersaturation dependence of the dif- 
fusion is an even 2nd order polynomial (aip 2 — b)/(a — b) 
the stable and unstable profiles still approach a lin- 
ear profile which approaches the spinodal concentration 
tp s = \fajb exactly at the critical point 



Nr. 



1 



1 



(12) 



where the stable and the unstable profile merge and dis- 
appear in a saddle node bifurcation. The representation 
of the stationary profiles via flow diagrams provides us 
with more insight into the structure of the solutions, and 
how they might change upon changes of the transport 
process. 



C. Flow diagram 

The different types of stationary supersaturation pro- 
files ip(z) can conveniently be summarized by observing 
that in Eq. ((8]) the control parameter N = £A 2 /D can 
be absorbed into the length scale, for instance by choos- 
ing z -> z ee z/y/N = {D/Cs 1/2 z/K Q. The equation 
((HI) has no free parameter in that case, such that each 
choice of initial conditions (y(l); <p'(l)) leads to a unique 
stationary solution of the diffusion equation. In other 
words, the solutions of Eq. © at different N only dif- 
fer by a rescaling of the length scale. Consequently, all 
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solutions can be represented as a flow diagram in the 
Qprf/VW) plane (Fig. [J). 

The flow is symmetric with respect to reflections at 
both coordinate axes. The symmetry with respect to re- 
flections at the horizontal axis reflects the free choice of 
the orientation of the z axis. The inversion symmetry 
of the plot is a consequence of the symmetry between 
the phases in the tp 4 model. In Eq. (JSj) these symme- 
tries amount to invariance under z — > — z and tp — > — tp. 
The symmetries are a consequence of the assumption that 
there are equal amounts of A and B in the mixture, and 
that the free energy is symmetric with respect to inter- 
changing A and B. App. lAldemonstrates that our findings 
do not substantially change in the non-symmetric case. 
There are only changes of unstable trajectories residing 
entirely in the spinodal regime where \tp\ < tp s . Up to this 
change the flow is structurally stable (cf. |30J| ) upon vari- 
ation of the functional form of the source term (as long as 
it remains continuous and strictly monotonous) and the 
dependence of the diffusion on the supersaturation (as 
long as one replaces (3p — l)/2 by a continuous convex 
function which takes the value unity at the equilibrium 
compositions <p = ±1 and negative values somewhere in 
between) . 

In order to identify the trajectories in Fig. [4] with the 
profiles shown in Fig. [2j we observe that all solutions 
shown in Fig. [2] have in common that zq — and tp(0) — 
1. They differ in the choice of the initial slope <p'(Q)/y/N. 
There is a stable solution with slope tp'(l) slightly smaller 
than y/ N/3, and an unstable solution where the slope is 
initially somewhat larger. For stronger driving the initial 
slopes of the solutions ever closer approach y/N/3, and 
the solutions remain linear for longer times. 

In the phase flow the profiles all correspond to trajec- 
tories which start on the line tp = 1 (rightmost dashed 
green line) and move inward with a negative initial slope. 
If > p'(0)/VN > -1/V3 one obtains a stable profile, 
which ends on the line <p' = 0. For tp'(0) < — 1/V3 one 
obtains an unstable profile, which terminates with infi- 
nite slope at the vertical green line tp = 1/V3, which 
indicates the spinodal ip s . For each choice of N there is 
exactly one stable trajectory and one unstable trajectory, 
which fulfill these boundary conditions for z = 1. In ad- 
dition there are symmetry-related solutions in the four 
other quadrants of the diagram. 

In addition to the supersaturation profiles, which start 
at if = 1 there are additional stationary states where the 
supersaturation always lies within the spinodal regime. 
More complicated profiles can be build by concatena- 
tion of pieces of profiles, provided that (3<p 2 — l)tp'(z) 
remains continuous in order to guarantee a continuous 
diffusion current. At least in principle, sharp changes of 
ip are admissible because they are the hallmark of phase 
boundaries. However, even for negligibly small M the 
forth order derivative term comes into play here, because 
it suppresses steep gradients and jumps in tp. Only the 
presence of this term makes it possible to discuss stable 
and metastable diffusion profiles connecting A rich and 
B rich states. 



IV. MODE EXPANSION 

The relaxation towards the stable profiles is a dynamic 
process. In order to also gain insight into the structure 
of the stationary profiles and their bifurcations, we resort 
to a mode expansion 



. . 2 ^-v . njz 



(13) 



where the prime at the sum indicates that the sum is 
to be taken over all odd integers j. Moreover, we re- 
quire (p-j — —tpj such that this ansatz automatically ful- 
fills the antisymmetry of the composition profile (which 
we had required to ensure mass conservation) and the 
boundary condition of the diffusion equation ((H)). In 
other words it shows an inversion symmetry at the phase 
boundary 2 = 0, and vanishing slope at the top and bot- 
tom z = ±1. 

Truncating the expansion after the first few leading 
order modes will provide us with a finite-dimensional dy- 
namical system mimicking the temporal dynamics of the 
system. 



A. Equations of motion 

Inserting Eq. (|13p into Eq. §5§ and using trigonometric 
relations in order to express products of sines in sines (cf. 
App. [B]) leads to 



N 



7T 2 f 7T 4 M 2 j 4 



8 



16 



k 2 



+ 3 X) ( kl + t ) f^m-k-i (i4) 



1. Single Mode 

Taking into account only the modes j = ±1 leads to 

3 , 



TT 2 7T 4 M 2 



(15) 



It admits the stationary solutions tp = and the nontriv- 
ial solution 



<y5 = 



7T 2 M 2 

6 



8N 
3^ 



1/2 



which ceases to exist at the critical point 



iV« = - 



TT 4 M 2 

16 



(16) 



(17) 



Already in this simplest conceivable framework it is not 
possible to have a stable diffusion profile when the pa- 
rameter N characterizing the strength of the tempera- 
ture drift becomes too large. However, in this setting 
the amplitude tp of the supersaturation profile vanishes 
upon approaching A^ c rather than that the stable profile 
vanishes in a saddle node bifurcation as suggested by the 
findings summarized in Fig. [3] 
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FIG. 5: (Color online) The panels show the amplitudes (a) (p\ and (b) (p$ of the steady-state solutions of the two mode 
approximations to the concentration profile Eq. (|18[) . Stable solutions are indicated by a solid red line, and the unstable ones 
by dashed green line. For comparison the dotted purple line in panel (a) also shows the single mode approximation Eq. (I16|l . 



2. Two Modes 



Taking into account the two leading order modes j = 
±1, ±3 leads to 



7T 2 7T 4 M 2 \ 3 , 3 2 „ , . 

<Pi = -[N- y+ 16 / 91 ~ 2 yi + 2 V ' lipa ~ Vltp5 ' 

. Ar 9tt 2 81tt 4 M 2 \ 9 o „„ 2 27 , ,„ nl , 

<^3 = -(AT- — + — ^ + (18b) 



For the case of narrow interfaces M = its steady-state 
solutions are plotted in Fig. [5] 

There still is a solution whose amplitude continuously 
disappears upon approaching the critical parameter value 
7V C (1) = tt/8 ~ 1.2337 found in Eq. (JTTJ). Close to the bi- 
furcations its amplitudes are very close to the one of the 
single mode approximation (purple dotted line). How- 
ever, this solution is now correctly identified as an un- 
stable solution where the concentration lies within the 
spinodal region all the time. It corresponds to one of the 
closed orbits in the center region of Fig. 
Moreover, the two mode model admits two qualitatively 
new types of solutions. For all N > there is a sta- 
ble solution (p x =Q and <p 3 = § U - - gf? f ) ^ 
which amounts to the single mode approximation of a 
system of height A/3 rather than A. 
In addition, for small N there are two (symmetry-related) 
pairs of non-trivial solutions, which disappear in saddlc- 

(2) 

node bifurcations at Nc ~ 0.5031, which is remarkably 
close to the exact value N c = (V3 — l) 2 ~ 0.5359 ob- 
tained in Eq. (JT2J) (cf. Fig. [2J). Beyond this value of N 
diffusion is too slow to remove the supersaturation from 
a domain of size A, and the system approaches the stable 
fixed point with vanishing amplitude (f>\. 



The saddle node bifurcation at N c and the subsequent 
decay towards a solution with domains of size A/3 has a 
simple physical interpretation. It amounts to nucleation 
of new domains close to the top and bottom of the system 
where supersaturation is largest, such that the domain 
size is reduced to one third of its previous value and dif- 
fusive relaxation can again keep track with the driving of 
the system. When relaxing the constraint of translational 
invariancc in horizontal direction, nucleation would of 
course give rise to a cloud of droplets rather than a full 
two-dimensional layer of material. A discussion of the 
related instabilities and feedback-mechanisms lie outside 
the scope of the present work. First steps to their under- 
standing have been suggested in [23| , and a more detailed 
analysis will be given in forthcoming work. At this point 
we conclude that the two-mode model faithfully describes 
the breakdown of stable diffusion by a saddle-node bifur- 
cation. Even though the values of the supersaturation 
are reproduced only in a crude approximation (dashed 
lines in Fig. [3]) its prediction of the critical point is off by 
less than 7%. Adding more modes in first place gives a 
more realistic description of the stable solutions remain- 
ing for N > N c and in the completely analogous saddle 
node bifurcations leading to breakdown of diffusion in 
these structures, when N becomes still larger. 



6 



March 15, 2008 — 12:46pm 



Phase Separation under Ultra-Slow Cooling: Onset of Nucleation 



V. DISCUSSION 

In this paper we have studied the diffusive transport in 
binary mixtures in response to ultra-slow cooling. Start- 
ing from the generic 4 model of binary fluid mixtures 
we showed 

• that stable diffusion is only possible up to a critical 
parameter value Eq. (fT2j) where the stable diffu- 
sion profiles merges with an unstable profile, and 
disappears in a saddle-node bifurcation. 

• This prediction of the critical value applies when- 
ever the concentration dependence of diffusivity is 
well-approximated by an even 2nd order polyno- 
mial. 

• In the limit of sharp interfaces the steady state pro- 
files may be represented in the form of flow dia- 
grams Figs. |4] and [71 Since the flow is topologically 
stable, the findings are not restricted to the spe- 
cific model but are robust to changes in the phase 
boundaries and concentration dependence of the 
diffusivity. 

• An expansion of the model in a Fourier series shows 
that the picture also applies when the forth order 
derivatives are included in the diffusion equation 
(|5j) that are needed to account for the finite width 
of the interface between coexisting phases. 

• When approaching the critical parameter value of 
the saddle-node bifurcation the concentration pro- 
file becomes unstable with respect to nucleation of 
droplets, which will break the horizontal transla- 
tion invariance assumed in the present study. 

The two-mode approximation of the concentration pro- 
files faithfully describes the approach to the bifurcation, 
and may hence serve as a starting point to further ex- 
plore the nucleation of droplets and its feedback on the 
transport of supersaturation. 

We conclude the present study by estimating the max- 
imal cooling rates where a mixture can still catch up with 
the change of temperature in a system of lateral extension 
A ~ 1 cm. For a mixture of methanol and hexane the 
equilibrium diffusion coefficient at the considered tem- 
peratures is of the order of D ~ 10~ 5 cm 2 /s. In view of 
Eqs. (U]) and (joa|) the system can only smoothly follow the 
change of temperature when the cooling rate the cooling 
rate does not exceed a critical value 



AT 
~dt 



N r . D 



&r<fo A 2 



(19) 



By comparison with the phase diagram for the mixtures 
(c/. Fig. 5 in [22|) one finds that for all experimentally 
accessible temperatures the cooling rate must not exceed 
a few K/h (Fig. [6]). For other binary mixtures these 
numbers will be of the same order, and for polymer so- 
lutions where the diffusion coefficient is still three orders 
of magnitude smaller the cooling rate must not exceed a 
few Kelvin per month. 




5 10 15 20 25 30 35 

T [°C] 



FIG. 6: (Color online) The solid red line gives an estimate (|19[) 
where the saddle node bifurcation leads to a breakdown of 
continuous diffusive demixing in the mixtures of hexane and 
methanol considered in detail in Ref. [22|. The label 'sta- 
ble' marks the range of heating rates v where diffusion keeps 
pace with the generation of supersaturation by the change of 
temperature. For larger heating rates (region marked as 'un- 
stable') supersaturation will become so large that droplets are 
nucleated. 



in view of the results summarized in Fig. [5] this study 
strongly suggest that complex time-dependent nonlinear 
behavior must be considered as the typical response of 
a binary liquid which is driven into the coexistence re- 
gion by a temperature ramp. Only when special precau- 
tions are taken, or when one approaches a region in the 
phase diagram where temperature changes do not induce 
changes of the composition diffusion can keep the system 
close to equilibrium. 
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APPENDIX A: ASYMMETRIC PHASE 
DIAGRAMS 

In (22| it is shown that binary mixtures with a sym- 
metric coexistence region are described by Eq. ([8|), even 
when the critical point does not correspond to equal vol- 
ume fractions of both phases. An asymmetric miscibility 
gap can be characterized by the temperature dependence 
of the mean value <fr(T) of the coexisting compositions. 
It gives rise to an additional source term £ = 0q 1 dt<j> in 
the diffusion equation, which breaks the ip — > — tp sym- 
metry. Absorbing £A 2 /Z? again into the length scale and 
concentrating on the case M <C 1 (i.e., S <C A) one 
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-1 
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the symmetric case Fig. [4] the heteroclinic connection be- 
tween the saddle points on the spinodal disappear. This 
leads to a rearrangement of the flow in the unstable re- 
gion — ip s < ip < ip s . On the other hand, there are only 
minor changes in the profiles outside this range, and the 
criterion |y'(0)| < 1 to obtain a stable diffusive profile 
remains valid to a very good accuracy. 



APPENDIX B: MODE EXPANSION FOR THE 
DIFFUSION EQUATION 



FIG. 7: (Color online) Flow diagram of Eq. (|A1|I for £/£ 
0.067. All other parameters are the same as in Fig. U 



When inserting the ansatz (|T3|) into the nonlinear dif- 
fusion equation (O only the nonlinear terms take some 
effort to evaluate. We observe that 



obtains 



3<^ 2 



d z ip 



(Al) 



The flow diagram for the stationary concentration pro- 
files for £/£ = 0.067 is shown in Fig. [7] In comparison to 



3^ 2 



d z ip 



3^ (d zV >) 2 + ^ 2 d 2 ^-\d 2 ^ 



and evaluate the terms one after the other. 

Observing that d z (p = ^ ■ j ifj cos we thus obtain 



2 6 v^'j • n j z n kz irlz 

3ip {d z <p) = - } y kl ifj tfk <pi sin — - cos — - cos — 



jki 



= ^- ^2 kl w fi 



an 



w(j + k + l)z Tr(j — k — l)z Tr(j + k — l)z n(j — k + l)z 
sin — — + sin — — + sin — — + sin — — 



Using the antisymmetry of tpj and shifting the summa- By an analogous calculation one finds 
tion index this yields 

3 2 n2 6 s-^'k 2 . irjz 

o fx n2 6 V-',, • n(j + k + l)z o z <p = - 2^ -7T Vj-k-lVk m sin— (Bib) 

3ip(d z ipy = - 2^ kl ifj tp k ipi sin ^ jki 2 

71 jki 

g , ■ Collecting the expansion coefficients of the Fourier series 

= — \ kl (fij-k-i tpk <Pi sin (Bla) then immediately leads to Eq. (|14p . 



jki 
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